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Abstract 

In this paper we continue the study of the three-layer Zamolodchikov model started in our 
previous works [18] and [19]. We analyse numerically the solutions to the Bethe ansatz equations 
obtained in [19]. We consider two regimes I and II which differ by the signs of the spherical 
sides (ai, 02, 03) — > (— aij — 02, —03)- We accept the two-line hypothesis for the regime I and the 
one-line hypothesis for the regime II. In the thermodynamic limit we derive integral equations 
for distribution densities and solve them exactly. We calculate the partition function for the 
three-layer Zamolodchikov model and check a compatibility of this result with the functional 
relations obtained in [19]. We also do some numeric checkings of our results. 
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1. Introduction 



This paper continues the investigation given in our previous papers [18], [19] of the three-layer 
Zamolodchikov model [1],[2]. In paper [9] Baxter considered the Zamolodchikov model with some 
special boundary conditions. Let us remind some basic ideas used by Baxter in [9]. To calculate 
the partition function Baxter used the symmetries of the Boltzmann weights of the Zamolodchikov 
model and some its special property, namely, the so-called " Body-Centered-Cube" (BCC) structure. 
The last property allowed to write down the Boltzmann weights up to some gauge and multiplication 
factors in the following form 

W(a\efg\bcd\h) = ]T ^^Km+km+k^k^) (L1) 

with all spin variables a . . . h taking only two possible values ±1 and <f) at b = —1 if a = b = —1, <p a ^ = 1 
otherwise. The sum in (1.1) is over the spin a = ±1 placed in the center of the cube corresponding 
to the weight W(a\efg\bcd\h) as it is shown in Fig. 1: 




e d 
Fig.l 



The coefficients Ki, K 2 , K 3 and K 4 were written in terms of parameters depending on the three 
angles &i, 9 2 , 3 and three spherical sides a±, a 2 , a 3 of the spherical triangle: 

v i = tanh2K i (1.2) 

where 

v x = -zT x T 2 , V2 = -iz^, v 3 = -z~ 1 T 1 T 2 , v 4 = iz' 1 ^ (1.3) 

and 

z = exp (ia 3 /2), T t = [tan (9 l /2)] 1/2 . (1.4) 

The BCC-form allowed to write down the transfer matrix T of the Zamolodchikov model as a product 
of two matrices 

T <t>,<t>' = XI ^W" Y <t>",<t>' (!- 5 ) 
<p" 

where 

Xw= II <f>«,c h e a{K3dh+K * ce \ Y r #= I] ^,a f e^ a9+K " bf) (1.6) 
cubes cubes 
and the summation in (1.5) is over central spins 0" between two successive layers of the lattice with 
the spins and <fi' respectively. Up to some simple multipliers defined by the formulae (6.7-6.8) and 
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(7.2-7.4) of [9] the Boltzmann weight W is equal to some other weight W which is invariant under 
the action of the group T of cube symmetries. This group has two generating elements r and p acting 
on the set of spins {a\ef g\bcd\h} (see, also, [17]) as follows 

r{a\efg\bcd\h} = {a\feg\cbd\h}, p{a\efg\bcd\h} = {g\cab\fhe\d}. (1.7) 

Under these transformations the spherical angles also change 

(91, 9 T 2 , 91) = (9 2 , 9 U 9 3 ), (91 91 9 P 3 ) = (tt - 9,, 9 3 , tt-9 2 ). (1.8) 

Let us go back to the BCC-form of the weight (1.1) and the expression (1.5) for the transfer matrix. 
Unfortunately, two matrices X and Y do not commute with each other. In order to handle this 
difficulty Baxter modified the boundary conditions for the central spins of Boltzmann weights so 
that each central spin a was replaced by the product of a pair of spins pp' as it is shown in Fig. 2 




e d e d 

Fig.2 



The matrices X and Y become 

X*#.= II ^,che^' iK3dh+K ' ce \ V,0'= II ^>/e^' (i ^ +X2b/) . (1.9) 
cubes cubes 

Now there exist some non-singular matrices P,Q,R,S such that 

X(v 3 ,v 4 ) = P(9 1 )R(v 3 ,v 4 )Q- 1 (9 1 ), Y( Vl ,v 2 ) = Q(9 1 )S(v 1 ,v 2 )p- 1 (9 1 ) (1.10) 



and the matrices R(v 3 ,v 4 ) and S(vi, v 2 ) are diagonal. As was mentioned above it leads to a mod- 
ification of boundary conditions for the initial spins a, namely, the product of all cr's along the 
front-to-back direction should be 1: 

aa a ... = pp p p p p ... = 1, (1.11) 

while the initial cr's were not constrained at all. If the cubic lattice is infinite in all three directions 
one can assume that the modification of the boundary conditions does not change the partition 
function. Further Baxter used two symmetries which do not affect front-to-back direction of the 
lattice. The first one is just the transformation r from (1.7). The second one, say r', is the reflection 
interchanging left and right 

T'{a\efg\bcd\h} = {f\dab\ghe\c} (1.12) 

with 

(9l',9 T 2 ',9l') = (n-9 1 ,9 2 ,n-9 3 ). (1.13) 
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Using these symmetries Baxter derived two functional relations for the largest diagonal elements of 
the matrices R and S in (1.10) and found its solution. It allowed him to come to the marvelous 
result for the partition function per site k (formula (8.14) and (9.7) of [9]) 3 : 

ln/c/2 = J(Ci) + J(C 2 ) + J(Cs) + J{Q (1-14) 

where 

J{Q = {2tiY 1 { (ln(2cosx) + xtaxix)dx (1.15) 

JO 

and 

e *i = e iCa = iv 2 , e iCi = -iv 3 , e KA = i/v 4 . (1.16) 

In paper [6] Baxter obtained the manifest result for the isotropic case 6i = a* = 7r/2 

« = 2 3 / 4 (2 1 / 2 - 1) e 2G/ " = 1.2480... (1.17) 

where G = 0.915965... is Catalan constant. In papers [7], [8] Baxter and Forrester checked this result 
numerically using the variational method. 

In [9] Baxter also obtained a generalization of the result (1.14-1.15) to the case of the lattice which 
is infinite in the vertical and left-to-right directions and has n layers in the front-to-back direction. 
The claim is that the formula (1.14) is still correct but the function J given by (1.15) has to be 
replaced by the function J n for which Baxter found the following formula: 

J„(C) = i»(C) + ^tan(| + ^) (1.18) 

with the function J n 

roo sinh 2xC 

J 4 — r_ h~ ( cosec h(2irx) — n (cosech(2irx/n))dx. (1-19) 

Baxter also checked that this result was valid for three particular cases, namely, n — 1 when the 
model was trivial, n = 2 when the model became the planar free-fermion model and for n = oo i.e. 
the solution (1.18-1.19) reduced to the formula (1.15) in the limit n — > oo. 

In this paper we would like to examine this result for the case n = 3 using the Bethe ansatz 
technique developed in our previous works [18], [19]. The formula (1.18) turns out to be incorrect 
for n — 3. We think that the main reason is that in the case when n = 3 one can not apply both of 
the symmetries r and t' to derive the functional relations and the resulting formula (1.14). As we 
show in section 2 the symmetry r' fails for n — 3. We expect that it also happens for any odd n > 1. 
Other reason is that the function J n should have more complicated analytic properties in comparison 
with that assumed by Baxter, namely, that J n (() is analytic in the strip |Re(£)| < ir/2. In section 9 
we discuss the analytic properties of our solution and check the validity of the inversion relation. 

One should note that the case of even n seems to be different. We have checked numerically that 
for n = 4, 6 both of the symmetries r and r' work. So we expect that they also work for any even n 
and the symmetry method by Baxter gives the right answer. 

In paper [14] Baxter and Bazhanov found a generalization of the Zamolodchikov model to the 
case of the arbitrary number of states iV and obtained the equivalence of the n-layer model with 



5 Here we have omitted the normalization function £ defined by formula (7.2) of [9] 
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modified boundary conditions and the s/(n)-chiral Potts model at q 2N = 1, [11], [12]. In particular, 
the modified three-layer Zamolodchikov model appeared to be equivalent to the s/(3)-chiral Potts 
model at q 2 = — 1. 

The main purpose of this paper is to calculate the partition function of the modified three-layer 
Zamolodchikov model or the s/(3)-chiral Potts model at q 2 = — 1 in the thermodynamic limit using 
the Bethe ansatz technique. We also perform extensive numerical calculations to support our results. 

The paper is organized as follows. In section 2 we discuss the symmetry properties of the modified 
three-layer Zamolodchikov model with respect to the r and r' transformations. In section 3 we es- 
tablish an exact correspondence between the three-layer Zamolodchikov model and the homogeneous 
s/(3)-chiral Potts model. In section 4 we discuss the Bethe ansatz equations corresponding to the 
largest eigenvalues of the transfer matrices. In section 5 we consider some special properties of the 
transfer matrix spectrum. In section 6 we introduce two-line hypothesis of solutions to the Bethe 
ansatz equations. We also derive and solve the integral equations for the distribution densities in 
the thermodynamic limit. In section 7 we calculate the partition function in thermodynamic limit 
using the solution for the distribution densities. A compatibility of the results from section 6 with 
the functional relations is established in section 8. In section 9 we consider the mechanism how the 
inversion relations for the partition function work. In section 10 we consider one-line regime. In the 
last section we discuss the results and outline some further developments. 

2. The r and r' transformations for the three- layer Boltz- 
mann weight 

As was explained in [14] the Boltzmann weights for the n-layer Baxter-Bazhanov model with N 
states and the s/(n)-chiral Potts model at q 2N = 1 can be connected with each other (see formulae 
(2.20) and (2.21) of [14]). In particular, for the three-layer Zamolodchikov model we have 

3 

W zam ({a}, {/3}, {7}, {6}) = n^k,7^ m |7m,am,AIA+i) (2.1) 

i=i 

where {a} = {aii, ol-i, CK3}, {f3} = fa, P3}, ■ ■ ■ Here we imply periodic boundary conditions in the 
front-to-back direction and each W % is given by (1.1). The Boltzmann weight (2.1) can be decomposed 
as follows 

W zam = W + + W. (2.2) 
where W + corresponds to the modified model described above with the constraint 

aa'a" = 1 

on the central spins of the weights and the weight W_ corresponds to a different choice of sign 

aa'a" = -1. 

Using formulas from the previous section we can write 

W ± ({a},{[3},{ 1 },{6}) = 
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EA A. A. p i^4(oia2MlA 1 2 + 0203M2A13 ± 0301/^3/^1) 

Yen ,±/ti/t3 Y012 ,M2Mi r"a3,/i3/i2 c 

M 

A , A ^ p ^3(A/32/tl/i2 + /3 2 /33/t2M3 ± A/3lM3/*l) 

rk A. As p -f^2(7l72/il/i2 + 7273M2M3 ± 737lM3Ml) 

^71 i/"l/"2 V72 ,/12M3 r73,±M3Ml C 

A c A r A r p -fsTl(<Sl<52/tl/i2 + <52<53/*2/i3 ± <53<5lM3Ml) 

V y fll,/tl/t2 l r()2,/t2/t3 l rfl3,±/'3Ml • v z "'-V 

Now let us consider the r-transformation from (1.7). Under this transformation {a} <-> {7} and 
#1 <-> #2 or i^i, i^ 2 , K 3 , K 4 — > i^i, — 27r/4, fT 3 , + 27r/4. Using the following relation 

A. A. A. p -jf (0102MIM2 + 0203M2/i3 +e 030lM3/ti) _ 

Y^ai,/ti/t2 l ra2,/t2/t3 Ya3,e M3/H c 

A> A> A> p ij(ena 2 + a 2 a :i +t a 3 ai - ^i/i 2 - M2/*3 - (2-e) /t 3 /ti - (l+2e)) (9 A\ 

t Ma 2 ,/ii/i2ra3,/t2M3 t rai,e/t3Mi c l^-^V 

with e = ±1 which can be directly verified we get the r-transformed weights W±: 

w^(K}, m, m, m) = ^^ w± a a }, m, { 7 }, {*» (2.5) 

where K}, {/T}, { 7 r }, {<T} = {7}, {a}, {5} and 

F({«}) = g*f ( a l a 2+0203+a3Cn) 

From (2.5) we see that W+ and W- transform in the same way. The ratio of the functions F is the 
gauge transformation, i.e. up to the same gauge transformations both weights W+ as well as W- are 
r-invariant. Of course, their sum W zam = W + + W- is also r-invariant. 

As for the second transformation r' given by (1.12) the sum W + + W_ is again invariant (up to 
some gauge) under r' . However, each term W± is not invariant under r' separately. It means that 
for n = 3 only the first r-symmetry is valid for the modified model but not r'. We have checked 
numerically that the same situation takes place for n — 5. We can expect that for any odd n > 1 
the symmetry t' is not valid. Therefore, for odd values of n > 1 one can not use both symmetries 
r and r' for a calculation of the partition function. However, as was mentioned in the Introduction 
both of the symmetries r and t' seem to work for any even n. 

3. The equivalence to the s/(3)-chiral Potts model 

As was shown in [14] the modified model with the Boltzmann weights W+ corresponds to the sl(3)- 
chiral Potts model at q 2 = — 1. In [18], [19] we used the expression for the "star" weight defined 
in [13]. It differs from the definition of [14] (see the formula (3.14) there) by a rotation of the star 
weight over the angle | round the front-to-back direction. However, for the modified model with 
n = 3 the weight W+ is not invariant under this rotation (as in the case of r'). 

Therefore, we need a direct connection of the Boltzmann weight for the modified model W + given 
by (2.3) with the definition (2.17) from [19] 

W°>» M% 6; P ,p>, g , q ') = ^ p ,(a,»)WM-»,»m p{ »,0) t (31) 

- Wp> q {8,H) 
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where the rapidity picture is defined by Fig. 3: 
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W qp (a,P) 



Fig. 3 



In (3.1) all indices are twofold, for example, a = (0:1,012) where 01,02 = {0, 1}. Let us denote 
these indices as Oj, fa etc. A correspondence with indices taking values ±1 is very simple: Oj = h(aii), 

A direct calculation shows that 

+R +R +R j ^-l^K^2)+h{S2&3)+h{S 1 8 2 )h{8 2 5 3 ) 
W+lW, {Pi, ill, W) = 2e 1 2 3 (_]_^(7i72)+^(7273)+ft(7i72)fc(7273) X 

^ star ({/i(oi0 2 ), /i(o 2 o 3 )}, {MAft), M&ft)}, {M7i72), ^(7273)}, {/i(5i5 2 ), q, q') 

(3.2) 

where the coefficient 2 is due to the fact that in (3.1) a summation is over the twofold index while 
in (2.3) indices are threefold. In addition, the parameters p,p', q, q' should be so that 

^ = 1 + v i = i-p'/q 6 ak 2= 1 + v 2_ i+p'/q' 



and 



1 -v 1 l+p'/q' l — v 2 1 - p'/q' ' 

e 4x 3 = 1 + ^3 = 1 + q'/p e 4x 4 = 1 + ^4 = 1 ^ 3 3 ^ 

i-v z 1 - <?7p' 1 - 1> 4 i-p/q 

q q' p p f a A \ 

vi = — 7, v 2 = — v 3 = — v 4 = - (3.4) 
p p q q 

with the parameters Vi given by the formula (1.3). Without loss of generality we can choose 

q=l, q' = -i/Tl p = v 4 , p' = -- (3.5) 

Vi 

in accordance with formulae (3.4) and with 7\ defined in (1.4). 



4. The partition function and Bethe ansatz equations 

Our main goal is to calculate the partition function Z for the modified Zamolodchikov model with 
three layers at least in the thermodynamic limit. If the lattice is finite and has N sites in the left- 
to-right direction and M sites in the vertical direction (in the front-to-back direction we have three 
sites) then 

Z = '£t + (K 1 ,K 2 ,K 3 ,K^ (4.1) 
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where t + (K 1 , K 2 , K 3 , K±)i is the J-th eigenvalue of the iV-site transfer matrix for the modified 
Zamolodchikov model with the weights W + given by (2.3) and the sum is taken over all eigenvalues. 
If the number M tends to infinity only the eigenvalue with maximal absolute value contributes (see 
[3]) and for the logarithm of the partition function per site k = Z X ^ NM ^ we have 

In K = _L i n t+ ( Ku k 2) k 3 , K J Q . (4.2) 

Now we have to connect the eigenvalues t+(Ki, K 2) K 3 , K A ) with the eigenvalues t(p; q, q'), t(p'; q, q') 
of the iV-site transfer matrices T(p\ q, q') and T(p'; q, q') of the s/(3)-chiral Potts model which were 
written in our previous work [19] as follows 

2 N 2 N 
t(p;q,q') = {p - g)N {p - ql)N s(p;q,q'), t(p;q,q') = {p _ q)N {p + ql)N s{p-q,q ! ) (4.3) 

where the functions s(p; q, q') and s(p; q, q') are polynomials in p of the degree k 

k k 

s(p;q,q') = a k (q,q') J| (p - p t ), s(p;q,q') = a k (q,q') ]J(p - ft) (4.4) 

and Pi, Pi are two sets of k numbers depending only on q, q'. Due to the correspondence (3.2) we have 
t+(K u K 2 , K 3 , K A ) = 2 N e 3N ^ +K * +K * +K Mp; Q, «W; 1, q')- (4-5) 

Hence, 

In («/2) = + K 1 + K 2 + K 3 + K 4 + ^ \nt(p; q, q') t(p'; q, q') . (4.6) 

where t(p; q, q')o and t(p'; q, q')o are the eigenvalues of the transfer matrices T(p; q, q') and T(p'; q, q') 
which give the maximal absolute value of the product t(p; q, q')ot(p'; q, q')o- As was discussed in [19] 
the degree of the polynomials k can be only 2N or 2iV — 1. In fact, we need to know only the product 
of unknown functions a n (q,q'),a n (q,q') which enter in (4.4), namely, 

a 2N (q,q') a 2N (q,q') = 4, a 2N _i{q ) q')a 2N _ 1 {q,q')=N(q' 2 -q 2 ) (4.7) 

to determine uniquely the eigenvalue t+(Ki, K 2 , K 3 , K±). The parameters pi and p { should satisfy 
the Bethe ansatz equations (5.4-5.5) of [19] 

and 



, , -Mf-'n^if (4.9) 



where 



f(p,x,q) = P —-^-. (4.10) 
p + q 



and u = e 2l7T / 3 is the root of unity of power three. 
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Below we shall consider two regimes I and II. Let us remind that the Boltzmann weights depend on 
spherical angles (#1, 9 2 , #3) and spherical sides (ai, a 2 , a 3 ) via the parameters Vi given by the formulae 
(1.3) and (1.4). Let us say that this parameterization corresponds to the regime I. 

The second regime II will correspond to the negating of spherical sides i.e. (01,02,03) — > 
(— Oi, — a 2 , —03) with the spherical angles {6i,6 2 ,6 3 ) being unchanged. In formulae for Vi only one 
substitution should be taken z — > z~ l i.e. new v\ are given by 

v [ = - Z ~ X T X T 2 , v' 2 = -iz- 1 T 2 /T 1 , v' 3 = -zT x T 2 , v' 4 = izT 2 /T ± . (4.11) 

However, we do not have a good geometric picture in regime II. 



5. The structure of the spectrum of the transfer matrix 

Let us introduce some notations we use below. There are three mutually commutative operators 
X 1 ,X 2 ,X 3 which also commute with both transfer matrices T and T 



where 



[X h T{p-q,q')\ 




[X h T{p-q,q')\ = 


: 


with 


Xi X{ 


<8> x 


i ® 




<8> x 




/0 


1 


°\ 




(0 


1 


0^ 




(0 








1\ 










1 




1 
















1 





Xi = 


1 








x 2 = 








1 


x 3 = 





1 










VO 1 





0) 




VO 





1 0) 




Vi 








0/ 



In addition, the transfer matrices T and T commute with shift operators P 2 and P 3 which act in 



left-to-right and front-to-back directions respectively. However, only the operator 

X = X 1 + X 2 + X 3 (5.3) 

commutes with both P 2 and P 3 . As a result all eigenvectors split into two big classes ip + and ip~ 

X^ + = i = 1,2,3, XiP + = 3ip + , Xip~ = -ip- (5.4) 

P 2 ^± = ftP^, n N = 1, p 2 = 0, N- (5.5) 
P^ ± = u p ^ ± ) co 3 = l, p 3 = -1,0,1 (5.6) 



It is interesting to note that all eigenvalues corresponding to the vectors ip + are in the sector 
with the degree k = 2N of the polynomials s(p) and s(p) and non-degenerate in general case. In 
contrary, the eigenvalues corresponding to ip~ are in the sector with k = 2N — 1 and have multiplicity 
3 in accordance with three eigenvalues of P 3 , p 3 — — 1, 0, 1. So let us use the following signature for 
eigenstates with fixed momentums p 2 and p 3 : (p2',Ps) + for ip + an d {p2\Pz)~ f° r ^ '■ 

Below we give a complete spectrum of the transfer matrix T(p; q, q')T(p'; q, q') for the cases 
= 2, 3, 4 and largest eigenvalues in sectors (0; 0) + and (0; *)~ (hereafter * stands for the values of 
Pz — — 1, 0, 1) for the cases iV = 5, 6, 7, 8, 9 with a special choice of the spectral parameters 

a = l, q' = -i, p = ef, p' = e - in/4 . (5.7) 
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It simplifies drastically the structure of the spectrum and corresponds to the symmetric point due 
to (3.5) 

7T 7T 7T 

*i=2» ^=2' 03 =2- (5 - 8) 
Numerical calculations show that such a choice does not change a multiplicity of the largest eigenvalue 
but it leads to some additional degeneracy in low eigenvalues. In general case the number of different 
eigenvalues grows too rapidly with N. 

In tables 1,2,3 we present the absolute eigenvalues A 

A = \t(p;q,q')t(p';q,q')\ (5.9) 

of the transfer matrix T(p; q, q')T(p'] q, q') for the cases N = 2, 3, 4 respectively. The upper subscript 
of eigenvalues shows its multiplicity. 





Ai 


A 2 


(0;0)+ 


16 a 




(0;±1)+ 


80 




(0;*)- 


16(3 + 2^) « 93.2548 


16(3 - 2^) « 2.7452 


(i;*)- 


16 2 





Table 1. The case N = 2. 





Ai 


A 2 


A 3 


A 4 


(0;0)+ 


2137.2 


38.81 


32 2 




(0;±1) + 


96 








(1;0)+ (2;0)+ 


96 








(1;±1)+ (2;±1)+ 


317.62 


250.91 






(0;*)- 


288 2 


119.43 


96 2 


8.574 


(1;*)- (2;*)" 


1120 


221.267 2 


13.884 2 





Table 2. The case N = 3. 



(P2;ps) e 


Ai 


A 2 


A 3 


A 4 


A 5 


A 6 


A 7 


A 8 


(0;0)+ 


5326.8 2 


1840.5 


1381.9 


256 


105.5 


49.21 2 






(0;±1)+ 


4920.0 


2720.5 2 


1456.2 


852.5 2 


135.8 








(l;0) + ,(3;0)+ 


3002.8 2 


164.77 2 














(2;0) + 


9741.8 2 


256 2 


242.2 2 












(1;±1)+ (3;±1)+ 


1858.2 


1313.8 


522.86 


229.6 










(2;±1)+ 


6269.5 


627.07 4 


381.53 












(0;*)- 


13548.7 


4518.5 2 


2304 2 


1062.1 2 


1024 3 


256 3 


127.6 2 


19.4 


(1;*)- (3;*)" 


7952.6 2 


3694.9 2 


1386.3 2 


603.3 2 


582.7 2 


467.1 2 


287.1 2 


28.2 2 


(2;*)" 


2458.9 4 


1536 2 


1511.2 2 


338.2 4 


271. 9 2 


40.8 2 







Table 3. The case N = 4. 
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In spite of the fact that there are complex eigenvalues in the spectrum, all eigenvalues in the sector 
(0; 0) + are real and positive. Bold numbers correspond to the eigenvalues with largest absolute value. 

In table 4 we give dimensions m and numeric values A = 7^ log A i only for the sectors (0;0) + 
and (0; *)~ with N — 2, ... ,9. 







N = 2 


N = 3 


N = 4 


N = 5 


N = 6 


N = 7 


N = 8 


N = 9 


(0;0)+ 


m 


2 


4 


8 


18 


66 


196 


696 


2440 


A 


0.4621 


0.8519 


0.7150 


0.7345 


0.8263 


0.7771 


0.7827 


0.8211 


(0;*)- 


m 


2 


6 


16 


52 


172 


586 


2048 


7286 


A 


0.7559 


0.6292 


0.7928 


0.8048 


0.7684 


0.8092 


0.8119 





Table 4. The largest eigenvalues in the sectors (0; 0) + and (0; *) . 



We have obtained these numbers by a direct diagonalization of the transfer matrix in correspond- 
ing subspaces with fixed values P2,P3,£ = ±1- The only missing number is for iV = 9 in the sector 
(0; *)~. The size of the transfer matrix becomes too large in this case. 

From these results we see that the largest eigenvalue belongs to the sector (0; 0) + with multiplicity 
1 for N = 31 and to the sector (0; *)~ with multiplicity 3 when iV is not divisible by 3. Of course, 
this degeneracy will not affect the value of A when the size of the transfer matrix goes to infinity. 
Moreover, the largest eigenvalues in the sectors (0; 0) + and (0; *)~ are strictly positive. Also we can 
see that when A" is divisible by 3 the value of A decreases and when A" is not divisible by 3 the value 
of A increases. One can assume that in the thermodynamic limit these two sequences tend to the 
same number. 

Numerical calculations show that for a different (close to the symmetric point) choice of the 
spectral parameters the picture remains essentially the same. We should note that when A" is not 
multiple of 3 this situation can change for a different choice of the spectral parameters. However, 
when N = 31 the largest eigenvalue is always in the sector (0, 0) + . Let us restrict our attention only 
to this case. As was mentioned above for the states in the class ip + the degree of polynomials s{p) 
and s{p) is k = 2N. 

Let us introduce the following parameterization for the parameters p iy p i satisfying (4.8-4.9) 

g-ivr/4 g«vr/4 

Pi = -Si^— , p i = Si^—, i = l,..,2N, (5.10) 
1 i 1 i 

Numerical analysis shows that Sj and are complex conjugate for all %. 

Moreover, for the states in the sector (0, 0) + 2A^ parameters s$ can be divided into two sets, 
namely, {sj, s i+ ^}, i — 1, N in such a way that 

SiS i+N = 1, i = l, ...,N. (5.11) 

It is not difficult to see that these equations reduce the system of equations (4.8-4.9). 
Let us write Sj and s« as 

Si = e Xi+iVi , s t = e Xi ~ iVi , % = 1, . . . , N, (5.12) 
where Xi and are two sets of real numbers. Let ( be a real parameter such that 

Ti = e- c . (5.13) 
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Then Bethe ansatz equations (4.8-4.9) can be rewritten in the following form 



cosh (xi + iyt ± i^) + cosh(C - i^) 
L cosh (xi + iyi ± i^) + cosh (C - i^) J 



. 7T s 



at cosh (xi + iyi ± i—) — cosh (xj — iyj) 

--n— — — f — - — - (5-i4) 



j=i cosh (xj + iy, ± «-) + cosh (xj - iyj) 



One has to note that the equations (5.14) should be valid for both signs simultaneously. Of course, 
its complex conjugations should be valid as well. So, in fact we have four sets of N equations. We 
can easily see that the transformation (xj,yj) — > (— Xj, — y$) is the symmetry of (5.14). 



6. The two-line hypothesis and the thermodynamic limit 

Now let us consider the structure of the solutions to the Bethe ansatz equations (5.14) corresponding 
to the largest eigenvalues in regime I when N = 31. In Fig. 4 we show the numerical data of the 
solution to the equations (5.14) for N = 31 and / = 1, . . . , 5 



y 



11.*-- 

12 71 

7T 

12. 


X 


11 -n-- 

12^ 

• • 


N = 3 

TT 

.12 X 


12" 

TT 

12. 


N = 6 


s •••••• 

• • • 


-ri * iV = 9 

12 " 

7T 

.12 . X 


• • • • 


iivr • iV = 12 

7T 

.12 X 




JV = 15 



Fig. 4. The structure of Bethe ansatz solutions for (xj,yj) in regime I 



All roots (xi,yi), i = 1, . . . , N split into two sets: 

lbr 2iV 7T 27V . . 

« = !,•• -,-3-; 2/» » ^2' l = ~ + 1 '---' iV - ( 6 - 1 ) 

For finite iV these two sets of roots form two curves in the complex plane. However, for large N 
these curves tend to two straight lines: y = ^ and y = -. Therefore, we expect that in the 
thermodynamic limit N — > 00 all solutions distribute on two lines with the imaginary parts ll7r/12 
and 7r/12. 

According to a standard method (see, for example, [3]) let us assume that in the limit iV — > 00 
the real parts X; L are distributed along the whole real axis with two different densities p_ and p + 
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for the first and the second lines respectively. Let Np_(k)dk and N p + (k)dk be the number of the 
solutions Xi with the imaginary part yi = Hit /VI and y^ = it j 12 in the interval k,k + dk. 

To write down the equations on the densities p + and p_ we need the Bethe equations (5.14) in 
the logarithmic form. Here one should control the phases coming from logarithms. Taking logarithm 
of (5.14) and making appropriate insertions of signs into the arguments of logarithms we can achieve 
the equidistant phase distribution for both lines. Then substituting the finite sums by integrals with 
the corresponding densities we arrive after some algebra at the following pair of the integral equations 



In 



s i n h(^-HL)sinh(^-HL) 



sinh(^ - f 



dk 



sinh(^ + ' 



+ ITT 



+ 



J°° dk'p+(k') lnj-tanh 



tanh 



fc-C 

2 1 12 

k — k! it 



k-k 1 



. TT 

12 



r2 

3 ~~ 

tanh 
k 



dk'p-(k') 



k + k' 



tanh 



2 
k' 



— 1 


TT ■ 


12 


TT- 




*6- 


}- 



(6.2) 



In 



cosh(^ + g) cosh( *=£ + f ) 



+ 



cosh (i±C - g) C osh(^ + f ) 
dk' p + {k r ) ln< — tanh 

J — oo I 

J dk'p_(k r ) lnj-tanh 



— ITT 



1 

- - 2 

3 



k — k' it 
+ , 6 



2 

k-k' 



tanh 



dk'p + (k') 



k + k' tt- 
h % — 

2 12 



TT 

+ I — 

12 



tan 



k + k' 



TT 

- ! 6 



= 



with normalization conditions 



dkp + (k) 



1 

3 : 



dkp_(k) 



2 
3' 



(6.3) 



(6.4) 



Then differentiating both equations by and making the Fourier transform we obtain the solution 
to these equations 

P±{k) = 2 cosh [2^- k)] ±1 (6 - 5) 

Substituting these formulae into the initial integral equations (6.2-6.3) one can check that they are 
really satisfied. 



In Fig. 5, 6 we plot the value 



as a function of Xi with 



llTT 



~ ^ for iV = 90 and the 



N x i+l -Xi " ^ -"i y« — i 2 5 Vi 

functions p~(x), p+(x) at C ~ —0.58666. From these figures one can see that discrete distributions 
are very close to its continuous limits. 
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N x i+1 - Xi 
N = 90, 1 < i < 59 



-2.5 -2 -1.5 -1 -0.5 



0.5 



x 



Fig. 5. Discrete distribution of x, with Hi = ^rf and the function p_(x 




iV Xj+i - Xj 
N = 90, 60 < i < 89 



x 



-2.5 -2 -1.5 -1 -0.5 0.5 
Fig. 6. Discrete distribution of x« with Hi = J2 ano - * ne function p+(x). 



7. Calculation of the partition function per site n 

In this section we calculate the partition function k given by (4.6) in the thermodynamic limit for 
the case of regime I. First of all, using the formulae (4.3,4.4,4.6,4.7) we obtain for the finite N 

In («/2) = - < ^+K 1 +K 2 +K 3 +K 4 + l ^-± In [(p + g )(p + </)(?' - + 9 , )] + ^i + ^2 (7.1) 
where 

n JV i AT i JV i AT 

5i = ^E ln (p-^) + ^^ ln ^-^' ^^^Ein^-pJ + ^E^^'-p^)- (7-2) 

1=1 1=1 1=1 j=l 
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Now let us consider the thermodynamic limit and replace all sums in (7.2) by integrals with 
corresponding densities. As above we can face the problem of sign choice in arguments of logarithms 
in integrals. One can choose these signs comparing each integral with a numerical value of the 
corresponding sum for finite N. Proceeding this way we obtain 

/OO fOO 
dkp-(k) In (p + e C-W4 e k+ilHr/12) + / dkp+ (ty l n (p + e C—/4 e fc+Wl2) + 
-00 J—oo 

/oo roo 
dkp.{k) ln (p + e C-W4 e -fc-ui7r/i2) + / ^(fc) ln (p + e C-W4 e -fc-Wi2^ ( 7 3 ) 

-oo J—oo 

/OO /'OO 
dA;p_(A;) ln (p' - e C+W4 e *-iiWi2) + / ^p + (A;) ln (-p' + ef+^e*^/ 1 ^ 
-oo J — OO 

/oo roo 
dkp_{k) ln (p' - e C+-/4 e -fc+,iWi2) + / ^p + (fc) ln {-p' + e C+W4 e -fe+Wi2)_ ( 7 .4) 
-oo J—oo 

Now using the solution for the densities p + and p_ (6.5), making some shifts in the integrals and 
taking into account (3.5) we come to the following formulae for Si and s 2 

s ^ r dk ln ^ 4 ^ ~ v 3~ lek+iw/6 )} ^ f°° dk ln K (! + V3 _1 e fe+i7r / 3 )] 
Sl 7r 7-oo 2 cosh 2A; — 1 7r J-oo 2 cosh 2k + 1 

VS r°° In \- e -' m ^ (1 - v 4 e k+i7T / 6 )] V3 r 00 In [e~^/ 3 (1 + v A e k+i7T ' 3 )} 

~ J-oo 2cosh2A;- 1 + ~ J-oo 2cosh2A;+ 1 ' ^ ' ' 



^ \n[^-e- i ^{l-v 2 - 1 e k+i ^)\ ^ ln e~^ 3 (1 + v 2 ~ V + -/ 3 )] 



s 2 =^r d k^i -— 

7r J-oo 2cosh2A; — 1 n J- 



dk- 

oo 2 cosh 2k + 1 



y/3 r° M-^l- V^ 6 )] ln [1 (1 + Vl e fc+ ^ 3 )] 

X£ / rfA: ^ — - — + rfA;^ — . (7.6) 

7T J-oo 2cosh2A;-l tt ./-oo 2cosh2A; + l 

Using formulas (1.2) and (3.5) we rewrite (7.1) after some algebra as 

ln («/2) = F( Vl ) + F^ 1 ) + F^ 1 ) + F(v 4 ) (7.7) 

for 

< 9i < tt, < at < n (7.8) 

where 

7TT 1 1 1 

= -- - - In (1 + „) - - ln (1 - „) + -I(v) (7.9) 

and 

J(u) = I + (ve l7T/3 ) +/_(-we" /6 ) (7.10) 

with 

\/3 r 00 rfA; 



J ± (z) = ^ r °* ln(l + ze fc ) (7.11) 
±v ; tt i-oo 2cosh2A;± 1 v ; v 
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Taking into account that v 2 = v%, v 3 = v\ and 

I{v~ l ) = f( v ) + \niv~ 1 ) + i(f)(v) 

where 



I( v ) = I + (ve- lw/3 ) + 14-ve-™/ 6 ) 



and 



An 2l\ 71 

11 < argt> < — 

9 3 & 6 

8n n 

— if - < argv < n 
9 6 

n An 

-- it n < argv < — 



(7.12) 
(7.13) 

(7.14) 



one can check that the pure imaginary term — ^ in (7.9) compensates the imaginary part in the 
whole expression (7.7) for all values of parameters in the domain (7.8). So the final expression (7.7) 
is real as it should be. We can also note that the form of (7.7) is similar to the Baxter's formula (1.14) 
for n = 3 with the function J 3 instead of J defined by (1.18) and (1.19). In terms of the variables 
Vi the formula (1.14) would be ®(vi) + ^(-v^ 1 ) + ^(v^ 1 ) + $(-f 4 ) with <&(-ie~ ix ) = J 3 (x). A 
difference with (7.7) is the presence of the additional signs before v 2 l and t> 4 . Moreover, the function 
F given by (7.9) is different in comparison with the function $. 

Now using some integral formulae from [4] and taking into account the cuts of the integrals (7.11) 
we can obtain the following expression for the combination (7.10) 

7T 71 

< arg z < — 

6 & 6 

7T 7T 

g <arg.<- 

tt 2n 
- < aigz < — 

2n 5n 



I(z) 



-j(o) 


(z) 








7(0) 


(z) 


-±ln(l 


-uoz 2 ) 




7(0) 


(z) 


-±ln(l 


- ujz 2 ) + ln(l - 


z 2 ) 


< J<°> 


(z) 


-±ln(l 


- ujz 2 ) + ln(l - 


uj- x z 2 ) 


7(0) 


(z) 


-±ln(l 


- uoz 2 ) + \ In (1 


-UJ^Z 2 ) 


7(0) 


(z) 


-ln(l - 


^ 2 ) + |ln(l- 


uj- x z 2 ) 


7(0) 


(z) 


- In (1 - 


c^ 2 ) + ±ln(l- 


uj- x z 2 ) +ln(l 



— < wgz < — 
3 6 

5n 7n 
— < arg z < — 
6 6 
7n 3n 

— < argz < — 

3lT 117T 

y <arg,< — 



(7.15) 



Here I^(z) corresponds to the principal branch of the integrals in (7.10) 

J<°>(z) = In (1 + uo' 1 z) + \ In (1 - u^z)— In (1 - uuz)+^-[2U 2 (1-z 2 )-U 2 (1-ujz 2 )-U 2 (1-uj- 1 z 2 )} 

2 2 An 



for — 7r/6 < arg z < n/6 
where uj = e t2w ^ 3 and Li 2 (z) is Euler dilogarithm 

In (1 -t) 



00 n rz 

Li 2 (^) = E £ 2 =" / 



(7.16) 



(7.17) 
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Now we can compare the asymptotics of ln(/c/2) with our result. Namely, for the parameters 
0i = 2.541,02 = 1.1,03 = 1-3 we get from the formula (7.7) approximately ln(/t/2) = 0.2224078.... 
Fig. 7 shows a dependence of ln(/c/2) on N for this choice of the spectral parameters. 



0.242 

0.232 

0.2226 
0.218 



ln(re/2) 



iV = 3A; 



3 6 9 
Fig. 7. 



30 60 
Logarithm of the partition function in regime I 



— 

N 



From Fig. 7 we can also see that for N > 30 the two digits in ln(«;/2) ~ 0.22 are fixed. It is in 
an agreement with its limiting value 0.2224078... 

In the end of this section let us consider an isotropic case of the Zamolodchikov model when all 
spherical angles and sides 0j = a% = ir/2. Then v± — — — e* 71 "/ 4 and V4 — I/V2 — e* 71 "/ 4 . Using the 
formula (7.15) one can get 

I(e**») + I(-e^) = g + 1 In (VS + 1) + ^ - | (7.18) 

where as in (1.17), G is Catalan constant. Substituting it into the formula (7.7) we obtain the answer 
for the logarithm of the partition function in the isotropic case: 

In («/2) = — + - i n (^3 + 1) _ ^1 ~ 0.354760685. (7.19) 
9n 9 18 

8. The functional relations in the thermodynamic limit 

In the previous sections we have obtained the integral equations (6.2-6.3) using the two- line hypothesis 
resulting from our numerical analysis of the solutions to the Bethe ansatz equations. Here we would 
like to discuss a compatibility of the results obtained above with the functional relations (4.1-4.4) of 
[19] for the transfer matrix eigenvalues t(p; q, q') and t(p; q, q') or (4.9-4.11) of [19] for s(p; q, q') and 
s(p; q, q'). Let us remind the reader that the functional relations for s(p; q, q') and s(p; q, q') look as 
follows 

s(p) s(p) s(-p) s(-ujp) = 
^0 s(p) s(-Ldp) + Af s(-p) s(-ujp) + A^ s(p) s{ujp) + A^ s(-p) s(ujp) 

(8.1) 

and 

s(— up) s(—p) s(p) s(p) = 
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\'q s(-ujp) s(p) + A'f s(-ujp) s(-p) + X' 2 s(u;p) s(p) + X'z s(ujp) s(-p) 

(8.2) 

where for simplicity a dependence on q, q' has been omitted and 

A = (p +uq) (p +uj- l q) (P + q') (P - q') 

Ai = (p + to q') (p + u' 1 q) (p + q) (p - q) 
A 2 = (p ~ q) (P + q) (p - uq') (p - q') 

A3 = (p ~ q') (p + ^~Y) (p - uq) (p - q) (8.3) 

with A^ being obtained from Aj by the substitution q — > —5. The relations (8.1) and (8.2) played a 
key role for a derivation of the Bethe ansatz equations (4.8) and (4.9). 
Denote _ 

u{p; q,q)= hm — In -, u(p; q,q) = lim — In ^— -. 8.4 

n^oo N a 2N (q,q') n~*<x> N a 2 N\q,q ) 

Using (4.7) and the fact that the relation (8.1) depends only on the product of s(p) and s(p) we note 
that a,2N(q,q') and a2N(q,q') disappear from (8.1). 

Now let us start with the expressions (7.3) and (7.4) for si and s 2 - In fact, si = u(p) and 
s 2 = u(p') in (8.4), but now we consider them as functions of arbitrary arguments (we fixed and 
omitted a dependence on q, q'). 

It is convenient to extract lnp terms from both expressions (7.3), (7.4). Due to the fact that 
J^° 00 dkp + (k) + J^ y 00 dkp^(k) = 1, it will give 2 hip for the u(p) and for u(p). Then taking into 
account the definitions (7.11) up to 2 hip we get 

u(p) ~ J + (-e-^ /6 ) + J + (-e- J7r/3 ) + J_(--e-^ 3 ) + 7_(--e- i7r / 6 ) (8.5) 
p p p p 

uip) ~ J+(--e l7r / 6 ) + J + (-I e -/3) + /_(£ e <T/3) + j_(I e -/6) ( 86 ) 
p p p p 

where for simplicity we denoted z — 1/T X 2 = e 2 ^ being not confused with another z defined in (1.4). 

Now we can consider logarithms divided by iV of the lhs of the functional relations (8.1-8.2) using 
the formulae (8.5) and (8.6). The same thing can be done for each of the four terms in the rhs of the 
relations (8.1-8.2). After this we can study a dominance of different terms in the rhs of (8.1-8.2) in 
the complex plane of the variable p. In Fig. 8 we show a behaviour of the logarithms of the absolute 
values for each of the four terms in the rhs of (8.1) as a functions of argp at iV = 60 dividing the 
result by N. 
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Fig. 8. A behaviour of logarithms of four terms in the rhs of the functional 
relation (8.1) within a region < argp < 2ir in regime I 

From this picture we can divide the complex plane into the four domains where the different 
terms dominate. Namely, 

I the first term dominates in the domain I: 

— 7r/6 < argp < 7r/6 and 5n/Q < aigp < 7tt/6 (8.7) 

II the second term dominates in the domain II: 

7r/3 < aigp < 27r/3 and 4tc/3 < aigp < 5n/3 (8.8) 

III the third term dominates in the domain III: 

27r/3 < argp < 5n/Q and 57r/3 < argp < ll7r/6 (8.9) 

IV the fourth term dominates in the domain IV: 

tt/6 < argp < tt/3 and 7tt/6 < argp < 4vr/3. (8.10) 
In Fig. 9 these four domains are shown: 




Fig. 9. 
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So, in the thermodynamic limit for these four domains (8.7-8.10) we should have four different 
relations. Namely, for the region (8.7) the first term in the rhs of (8.1) dominates and we should 
have the following relation 

u(p) + u(—p) = In [(p + uq){p + u~ 1 q)(p + q')(p - q')}. (8-H) 

Extracting the terms like hip as was described above and using (3.5) we get for the lhs of (8.11) 

j /£ e -Mr/6N + j (I e -<T/3) + /_(_£ e -W3) + J_(_I e -W6) 

p p p p 

/ + (£ e -/6) + / (I e W3) + 7 _(_£ e W3) + 7_(_I e <T/6). ( 8 . 12 ) 

p p p p 

Now using the following relations 

I + {e in/2 w) + I_(w) = In (1 + e^ /6 w) for - tt < argw < tt/2 (8.13) 

I + (e t7r/2 w) + /_(«;) = In (1 - e" m/6 w) + In (1 + e~ t7r/6 w) - In (1 - e t7r/6 w) 

for 7r/2<argw<7r (8.14) 

I + (e- tn/2 w) + 7_H = In (1 + e- m/6 u>) for - tt/2 < argw < tt (8.15) 

J+(e- l7r / 2 w) + J_H = In (1 + e lw/6 w) + In (1 - e ln/6 w) - In (1 - e-^'^w) 

for tt < argw < 3vr/2 (8.16) 

one can get the following result for the expression (8.12) combining the pairs of I + and 7_ in accor- 
dance with (8.13-8.16) 

In (1 + -e w/2 ) + In (1 + - e t27T/3 ) + In (1 - -e l7r/2 ) + In (1 + i e -^/ 3 ) = 
p p p p 

In (1 + q - ) + In (1 - ^) + In (1 + -) + In (1 + — ). (8.17) 
p p p p 

It is easy to see that the last expression is just the rhs of (8.11) after extracting logarithms of p. For 
the second domain II (8.8) the second term in the rhs of (8.1) dominates and we should have the 
following relation 

u(p) + u{p) = In [(p + uq')(p + u- l q'){p + q)(p - q)} (8.18) 

which can be checked in a similar way as (8.11) using relations (8.13-8.16). For the third domain III 
(8.9) the third term in the rhs of (8.1) dominates and we have to check a little bit more complicated 
relation 

u{p) + u(—p) + u(—ujp) — u(ujp) = In [(p — q)(p + uj~ l q)(p — coq')(p — q')}. (8.19) 
In the lhs of (8.19) we get 

j (£ e -*T/6) + 7 + (i e ~ m/3 ) + I(--e- lw/3 ) + I(--e~ tn/6 ) + 
p p p p 

7 (£ e -/6) + 7 (I e W3) + 7_(_£ e "r/3) + 7_(_I e W6) + 

p p p p 
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J + (-e^ /6 ) + /+(-) + /_(--) + I_(--e™'«)- 
p p p p 

I + (-- e ™/ 6 ) - /+(--) - /_(-) - I-(-e^ e ). (8.20) 
p p p p 

Now using (8.13-8.16) we combine all pairs of 1+ and I in (8.20) and then obtain 

In (1 - -) + In (1 - -e tn/3 ) + In (1 + -e^ /3 ) + In (1 - -e l7r/6 ) 
p p p p 

which is the rhs of the relation (8.19) after extracting the terms hip. For the domain IV (8.9) where 
the fourth term in the rhs of (8.1) dominates we have the similar to (8.19) relation 

u(p) + w(p) + u(-up) - u(up) =]n[(p- q') (p + w~V ) (p - uq) (p - q)] (8.21) 

which can be checked in the same way as the relation (8.19). 

The second functional relation (8.2) can be proved similarly. So we have checked that the func- 
tional relations (8.1) and (8.2) are compatible with our result obtained from the two-line hypothesis. 
Therefore, this fact can be considered as a justification of this hypothesis. 



9. Inversion relations 

In this section we want to check a compatibility of our solution (7.7) with the inversion relations. 
Considering the product of two transfer matrices T p T p i one can see that this product is decomposable 
in two cases. Namely, one can fix the parameters as follows p' = p and p' = —p. It is easy to see 
that the star weight defined by the formula (3.1) satisfies the following property: 

W star (a,p n ,a;p,p,q,q') = % q ,5^ W star (a, [3, (3, S;p, -p, q, q') = % q 5 a5 (9.1) 

where the inversion factor 

4(p3 + q 3 ) 

Due to the property (9.1) one has the following decompositions of the product of the transfer matrices 

T p T p = ^ q/ (I + T^), T p T_ p = ^ q (I + T^) (9.3) 

where I is the identity matrix of the dimension 4^ and (T^) is a transfer matrix built up from 
some weights (W^). Unfortunately, at the moment we can not establish rigorously in what 

domain of the parameter p, q, q' the terms and in the rhs of equations (9.3) can be neglected 
when iV — > oo. So we want to consider how the inversion relations work for the answer (7.7). 

Let us start with the expression for the eigenvalues t{p) and t{p) in the thermodynamic limit. 
Using the results of section 7 we get 

1 2m 1 2m 

lim -lnt(p) = ln2+j(^ 1 )+jM+a+ — , lim -lnt(p') = \n2+j(v 1 )+j(v^ 1 )-a-— (9.4) 

N—>oo iv O N-*oo iv O 

where p, p', q, q' now are given by (3.5) and 



m 



j(v) = ---\n(l + v) + I(v) (9.5) 
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with I given by (7.10). The function a = lim.N^ oc l/Nlog(a 2 N{q,q f )) with a 2 N entering to (4.4) is 
unknown but it cancels in the sum lnt(p) + lni(p'). It is evident, that the function F defined by 
(7.9) is 

F(v)= 1 -\n 1 T ±^+ 1 -j(v) (9.6) 
and the logarithm of the partition function given by (7.7) can be also represented as 



ln( K /2) = K 1 + K 2 + K 3 + K 4 + ~x(vi,v 2 , v 3 , v 4 ) 



where 



X(vi,v 2 ,v 3 ,v 4 ) =j{vi) + j(0 +J'(^3 1 ) 
It can be checked using the formulae (8.13-8.16) that 



(9.7) 
(9.8) 



j(v)+j(-v x ) 




57T TT 

if — — < argv < -- 
6 6 

.„ tt 2tt 
if - < argv < — 



(9.9) 



j(v)+j(v )=ln 



1 + v 3 



+ < 



2m tt tt 

-— if -- < argw < - 
3 6 6 

AlTT . 27T 

-— if —<axgv<ir 

2i7T . r 2tt 

— if -7r < argf < — — 



(9.10) 



Let us remind that in the domains I and II defined by (8.7) and (8.8) the inversion relations (8.11) 
and (8.18) are valid for s(p) and s(p) and therefore for t(p) and t(p). Another way to prove them 
is to use relations (9.9-9.10). Indeed, if p' — p then l/v 3 = q'/p,v 4 — p, v\ = —1/p, l/v 2 = p/q' and 
using (9.4) with the relations (9.9,9.10) we get 

lim -J- \nt(p) + lim -J- \nt(p) = In 4 + j{p) + j(-l/p) + j(p/q') + j(q'/p) = 



In 4 + < 



2m 
2iir 



In 



+ In 



p° + q 



+ 



2m 



(p + q'f ' 3 
p 3 + q' 3 Am 
(p + q'f ~ ~ 



In <& pql + < 



2iir p 3 + q' 3 2m 

— + n (p + g')3 " — 

7T 27T 

—2m for — < argp < — 
for - < argp < - 



- IT IT 

tor — < arg p < — 
3 By 2 

r it 2tt 
for - < argp < — 

. 2tt tt 
— 3~ < aTgP < ~ 3 



2tt tt 
and — — < argp < -- 



i.e. up to 2iir the inversion relation (8.18) is valid in the domain II. 



(9.11) 
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For the second case p' = —p we have l/v 3 = q' /p, V4 = p,v± = 1/p, l/v 2 = —p/q' and we have 



lim -1 lnt(p) + lim -1 lnt(-p) = In 4 + j(p) + J'(Vp) + j(P/V) + j(-?7p) 

N—>oo iv N—KX) iv 



In 



7T 



7T 



In 4 + < 



In $ p9 + < 



p 3 + 1 2i7r 2m 

^TTr- — + — for -6 <ar ^<6 

p 3 + 1 4?7r 2m 5n 

— for Y <arg P <vr = 

. p 3 + 1 2i7T 2i7T „ 57T 

l ln (^TIF + — "X for < argp < - y 

-2m for — < argp < 7r 
6 

7T 7T 5/T 

for — — < argp < — and — n < argp < — — 
6 6 3 



(9.12) 



i.e. again the inversion relation (8.11) is valid in the domain I up to 2iir. 

As we show below the inversion relation can also be written for the partition function. If 
2n/3 < a 3 < n then for the arguments of parameters Vi given by (1.3) one has 



2tt 

T 

2tt 



7T 



7T 



— 3- < argf^i) < --, -- < arg(-t; 1 ) < 



and 



: , < arg(v 3 l ) < ~, ~ < arg(-v 3 ) < 



~ < arg(t> 2 ) < 0, < arg(t; 2 v ) < ^ 
■- < arg^J 1 ) < 0, < arg(w 4 ) < - 



71 

3 

7T 

3 



(9.13) 



(9.14) 



So since (9.13) is valid one can use (9.9) with v — vi, v — —v 1 1 , v = v 3 1 or v = —v 3 , namely, 

(9.15) 



J{vi) +j(-v 1 ) = —, j{v 3 ) +j(-v 3 ) = — . 



Similarly, due to (9.14) the relation (9.10) is satisfied for v = v 2 , v = v 2 \ v = V4 , v = v 4 1 and one 
has 

p 3 + q 3 2i7i 
3"' 



j{v 2 ) +j(v 2 



In 



(9.16) 



Since, the logarithm of the partition function (9.7) is connected with the function j via the combi 



nation \ let us show the inversion relation for the function x- 

X(vi,v 2 , v 3 , v 4 ) + xi-vT 1 , v 2 \ -v 3 \ vl 1 ) = In • 



3 1 3 /3 1 /3 

1 . . — 1 . , -k _ ,„ p - <r _ ln p ■ q 



(p + <?) 3 (p' + q') 3 



(9.17) 



So in the rhs of (9.17) we have got the sum of logarithms of the function which is proportional to 
the product of the inversion factor (9.2). 
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Let us note that the function j(—ie~ lx ) resembles the Baxter's function ^(x) (1.19) because 
it satisfies the relations (9.9) and (9.10) which are similar to the relations (10.1, 10.2) of [9] 4 A 
difference is that the relations (9.9-9.10) are defined only in the restricted region and have additional 
pure imaginary terms in the rhs which are connected with the cuts of the function I entering to 
(9.5). However, as we have shown above the relations (9.9-9.10) are enough for the validity of the 
inversion relation (9.17). The function Js(x) is analytic in a wider strip |Re(x)| < ir/2. Of course, it 
also ensures the validity of the inversion relation (9.17). However, the applicability domain for the 
inversion relation ensured by the function J^{x) is wider. Namely, it works for the region < 03 < n 
while our result satisfies the inversion relation when 2n/3 < a 3 < n as was discussed above. 



10. The one-line regime 

Now let us consider the second regime II which differs from the regime I by a transformation — > — a*. 
As was discussed above the parameters V; L should be substituted by the parameters v[ defined by 
(4.11). 

The numerical analysis shows that the transformation — > — a« changes the ground state dras- 
tically in comparison with the regime I where the two-line solution corresponds to the dominating 
eigenvalues of the transfer matrices. Now the one-line solution corresponds to the dominating regime. 

Namely, for the one-line solution the imaginary parts defined by the formula (5.12) are close to 
77r/12 for alH = 1, . . . , N. As in the two-line case a precision of the approximation — In /12 m 
becomes better with a growth of N. So, we also can accept the one-line hypothesis that in the 
thermodynamic limit iV — > 00 the line corresponding to i/i = In/ 12 tends to be exact. It is interesting 
to note that now it is not necessary to consider only those N which are multiples of 3 as it was in 
the two-line case. The one-line solution works for any N. 

Now we would like to generalize the result of the sections 5 and 6 to the one-line case. Somehow 
this case is simpler in comparison with the two-line case because now we are dealing only with one 
line of the roots. Therefore, it is not necessary to divide the interval 1, . . . , N into two parts as we did 
it in the section 5. It is not very difficult to repeat here a strategy of the section 5, i.e. to introduce 
a distribution density p of the real parts in the thermodynamic limit and to obtain the integral 
equation for it. As a result we have the equation 



In 



cosh (k — in/A) — cosh (^ — in/12) 
cosh (k — in/A) — cosh (( — ibn/Yl) 

cosh (k — in/A) — cosh (k' + i5n/12 



+ in (1 - 2 [ k dk'p(k')) 



/oo 
dk'p(k') In 
-00 



cosh (k - in/ A) + cosh (k' + i5n/12) 



:io.D 



with the normalization condition 



/oo 
dkp{k) = 1. (10.2) 
-co 



As the integral equations from the section 5 the equation (10.1) can also be solved exactly. The 
result looks as follows: 

P{ > 2 cosh [2(C + &)] + ! 2cosh[2(C-fc)]-l' 1 ' 



4 The relations (10.1, 10.2) of [9] were written for the function J„ (see (1.18)). To get such equations for the function 
J n one should take into account an additional term 1/4 lntan(x/2 + 7r/4) in (10.5) of [9] 
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Actually, one can rewrite this result as 

p(k) = p + (-k) + p_{k) (10.4) 

where p± are the densities in the two-line case given by (6.5). It is quite amazing that the density in 
the one-line regime appears to be a superposition of the densities for the two-line case. 

Now we can calculate the logarithm of the partition function per site k which is still given by 
the formula (7.1) of section 6. For the one-line regime Si and s 2 which were introduced by (7.3) and 
(7.4) for the two-line regime are defined by 

/OO fOQ 
dkp(k) In (p + e C-i*/4 e *+i7*/12) + / dkp ^ ln (p + e C-^/4 e -fe-^/12) 
-OO J — OO 

/OO fOC 
dkp(k) ln (p' - e C+W4 e fc-i7 W /i2) + / rfA;p(A;) ln ^ _ gC+^g-fe+^Z^) ( 1Q 5 ) 

where p and p' are still given by formula (3.5) but for parameters V{ substituted by v[ defined above. 
Using the solution (10.4) and taking into account (3.5) we get s\ and s 2 in the thermodynamic limit 

ln[-4e-" /6 (l-t>3e* +i " /6 )] m In (1 + v',e k+ "' 3 )] 



7r J-oo 2cosh2A; — 1 7r J-oc 2cosh2A; + l 

7T J-oo 2 cosh 2fc - 1 7T J-oo 2 cosh 2A; + 1 



l n [_J_(l_^ e fe+-/6)] ln [_J_(l + ^ e fe+-/3)] 

S2 = ^ r ^ ^ + ^ r ^ ^ 

7r 7-oo 2 cosh 2k — 1 7r i-oo 2 cosh 2k + 1 

V3 [°° In [e~ l7r / 6 (1 - ^~V +OT / 6 )] [°° In [-e"^ 3 (1 + ^"V+^Z 3 )] 

H / 1 1 ~\~ / Cbki : : . ( lU.Y) 

7T J-oo 2 cosh 2k - 1 7r 7-oo 2 cosh 2A; + 1 

After some algebra we come to the following result for the logarithm of the partition function in 
the one-line case: 

ln («/2) = F'{v'- 1 ) + F'(v' 2 ) + F>(v> 3 ) + F'(v'f l ) (10.8) 

where the function F'(v) differs from the function F(v) defined by the formulae (7.9-7.11) only by 
the pure imaginary term: 

F'{v) = ^ - 1 ln (1 + v) - \ In (1 - v) + ±I(v). (10.9) 

One can check that, as in formula (7.7), the pure imaginary term compensates the imaginary part 
coming from other terms making the final answer (10.8) real. 

Let us note that proceeding as in the section 7 we can also check that the one-line solution is 
compatible with the functional relations in the thermodynamic limit. 
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11. Conclusion 



Let us summarize the basic results obtained above. Our main point here was to calculate the partition 
function of the Zamolodchikov model for n = 3 layers based on the functional relations and the Bethe 
ansatz equations obtained in our previous works [18], [19]. We have considered two regimes I and II. 
Analysing numerically the structure of the solutions to the Bethe ansatz equations corresponding 
to the largest eigenvalues we have come to the two-line hypothesis for the regime I and one-line 
hypothesis for the regime II. Using the standard technique we have got the integral equations for 
the distribution densities in the thermodynamic limit which are solvable exactly due to a difference 
property of the kernel of these integral equations. We have compared the results obtained from the 
integral equations with the discrete distributions for the finite size of the lattice. A similarity of the 
discrete distributions with the corresponding solutions of the integral equations pretending to be the 
continuous limits of the discrete distributions tells us about a plausibility of our conjectures. 

We have used these distribution densities for a calculation of the partition function. The answer 
can be expressed in terms of the dilogarithm function. It is in a good agreement with the data 
obtained from the numeric solutions to the Bethe ansatz equations up to N = 90. We have also 
checked a compatibility of this result with the functional relations. 

It would be very interesting to apply the above calculations to the three-layer Zamolodchikov 
model without a modification of the boundary conditions. We hope to do that in our further publi- 
cation. We also hope that proceeding this way we could get a better understanding of the problem 
discussed above and to study in thermodynamic limit the ground state of the hamiltonian derived 
by Baxter and Quispel in paper [10]. 

We also intend to develope the thermodynamic Bethe ansatz (TBA) technique for a study of the 
finite size corrections and possible conformal properties of the three-layer Zamolodchikov model. 

As a further step we would like to generilize our results obtained for the three-layer case to a 
generic case of arbitrary number of layers. 
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